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A  3-D  Plasma  Interaction  Modeling  System  (PIMS)  is  being  developed  to  predict  the 
interaction  of  electric  propulsion  plumes  with  surfaces.  The  system  is  designed  to  be 
flexible,  usable,  and  expandable,  allowing  users  to  define  and  mesh  surfaces  with  their 
choice  of  off-the-shelf  3-D  solid  modeling  packages.  These  surfaces  are  then  loaded  into 
PIMS,  which  performs  plasma  operations  based  on  user  commands.  Functional  PIMS 
modules  will  range  from  simple  (prescribed  plume  field)  to  complex  (full  PIC-DSMC) 
depending  on  the  user’s  request.  PIMS  will  compute  surface  interaction  parameters  such 
as  ion  flux,  ion  energy,  sputtering,  and  re-deposition.  Development  of  PIMS  to  this  date 
has  progressed  to  include  modules  that  a)  import  and  superimpose  prescribed  plume 
distributions,  and  b)  perform  ray  tracing  of  flux  from  point  sources.  This  paper  presents 
some  of  the  first  PIMS  results  -  sputtering  predictions  on  a  spacecraft  and  in  a  vacuum 
chamber  due  to  a  Hall-effect  thruster. 


Introduction 

Onboard  electric  propulsion  (EP)  thrusters,  which  use 
electric  power  to  generate  or  augment  thrust,  hold  the 
promise  of  greatly  increased  satellite  maneuverability, 
and  enabling  new  missions.  Many  types  of  EP 
thrusters  are  already  in  mature  states  of  development, 
and  many  can  achieve  specific  impulses  over  3000 
seconds.  This,  combined  with  growing  electric  power 
levels  onboard  new-generation  spacecraft,  is  pushing 
EP  rapidly  into  the  mainstream.1,2 

Several  EP  devices  are  currently  being  evaluated  for 
use  onboard  U.S.  commercial  and  military  spacecraft. 
One  of  the  most  promising  for  near-term  use  is  the 
Hall-effect  thruster  (HET).  Over  120  HETs  have 
flown  on  Russian  spacecraft,  where  typical  flight  units 
have  specific  impulses  around  1600  seconds  and 
efficiencies  near  50%.3  HETs  operate  by  generating  a 
stationary  xenon  plasma  inside  an  annular  channel. 
Strong  radial  magnetic  fields  are  applied  which 
impede  electron  motion,  but  allow  ions  to  accelerate 


axially  out  of  the  device  with  velocities  around  20 
km/s  (energies  of  around  300  eV). 

High-energy  HET  exhaust  ions  may  erode  (sputter) 
surfaces  on  which  they  impinge.  In  addition,  this 
sputtered  material  may  be  re-deposited  on  other 
spacecraft  surfaces.  These  issues,  and  others,  such  as 
electromagnetic  interference  and  spacecraft  charging, 
cause  some  concern  for  spacecraft  designers  who  want 
the  maneuverability  EP  offers  but  do  not  want 
increased  risk. 

Efforts  are  underway  to  quantify  some  of  the  risks 
associated  with  integration  of  EP  with  spacecraft, 
including  surface  erosion  and  re-deposition.  Work  has 
been  done  to  computationally  model  expansion  of 
HET  plumes.4  Additionally,  Gardner  et  al.  have 
developed  Environment  Work  Bench  (EWB),  a  code 
that  calculates  sputtering  of  spacecraft  surfaces  by 
superimposing  pre-computed  EP  plumes  onto 
spacecraft  geometries.5  However,  existing  codes  do 
not  self-consistently  calculate  the  plume  expansion 
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Fig.  1.  PIMS  data  structure.  Each  module  represents  a  cohesive  block  of  code  and  data  with  a  specific 
function. 


with  the  3-D  surface  geometry  in  a  usable,  flexible 
way. 

This  paper  describes  the  architecture  and  function  of  a 
new  software  package  named  Plasma  Interaction 
Modeling  System  (PIMS),  which  is  being  developed 
collaboratively  by  the  U.S.  Air  Force  Research 
Laboratory  and  the  Massachusetts  Institute  of 
Technology  with  the  goal  of  self-consistently  modeling 
plume  expansion  and  interactions  with  arbitrary  3-D 
surfaces.  Three  important  requirements  have  been 
placed  on  PIMS:  It  must  be  USABLE,  FLEXIBLE, 
and  EXPANDABLE. 

USABLE  means  a  typical  engineer  should  be  able  to 
set  up  and  run  a  typical  low-fidelity  case  in  less  than 
one  day  with  less  than  three  days  training. 

FLEXIBLE  means  PIMS  must  be  able  to  simulate  at 
least  three  important  cases:  a)  a  single  spacecraft,  b) 
multiple  spacecraft  in  formation,  and  c)  laboratory 
conditions  (e.g.  the  interior  of  a  vacuum  test  facility). 
Simulating  laboratory  conditions  is  very  important  for 
two  reasons.'  First,  since  there  is  very  little  on-orbit 
data  for  EP  thrusters,  ground-based  tests  must  be 
relied  upon  for  the  bulk  of  code  validation.  Second,  by 
modeling  the  laboratory  conditions,  PIMS  can  help 
engineers  interpret  lab  measurements. 


EXPANDABLE  means  PIMS  can  be  easily  expanded 
to  incorporate  new  plasma  simulation  algorithms,  new 
capabilities,  or  improved  efficiency. 

Approach 

PIMS  has  been  designed  as  a  collection  of  modules, 
each  with  a  specific  function  and  hierarchy.  Each 
module  contains  data  and  associated  code.  Modules 
may  be  categorized  into  three  levels,  as  shown  in  Fig. 
1. 

Level  0  modules  perform  functions  related  to  user- 
interaction.  Although  PIMS  is  fundamentally 
command-driven,  a  Graphical  User  Interface  (GUI) 
front  end  is  envisioned  for  the  future. 

Level  1  modules  are  the  primary  components  of  PIMS. 
They  perform  functions  related  to  propagating  plasmas 
on  the  volume  domain.  They  contain  algorithms,  such 
as  fluid,  PIC,  DSMC,  or  hybrids  thereof,  which 
perform  a  solution  subject  to  pre-set  boundary 
conditions.  Level  1  modules  are  uniform  in  that  they 
all  conform  to  a  specific  Interface  Control  Document 
(ICD)  -  they  have  specific  inputs,  outputs,  and 
resources  available  to  them. 


Table  1.  Sputter  yield  coefficients  for 
bombardment  by  singly  ionized  xenon. 


Material 

Coefficient  a 

Coefficient  b  (J'1) 

A1 

1.0 

1.9el6 

ITO 

0.1 

6.25el5 

Kapton 

0.05 

2.5el4 

AgT5 

1.0 

1.9el6 

Level  2  modules  set  boundary  conditions,  and  provide 
support  to  Level  1  modules.  They  act  as  a  toolbox  or 
collection  of  resources. 

The  purpose  of  the  modular  design  is  to  give  PIMS 
flexibility  and  expandability.  A  large  number  of  Level 
1  modules  are  desired  to  allow  flexibility  in  solving  a 
variety  of  different  problems.  The  ICD  is,  therefore, 
very  important,  because  it  describes  for  authors  of 
Level  1  modules  a)  what  inputs  and  boundary 
conditions  must  be  recognized,  b)  what  outputs  are 
expected,  and  c)  what  Level  2  resources  are  available. 
The  ICD  may  be  distributed  to  outside  groups  so  that 
PIMS  can  be  expanded  through  addition  of  new  Level 
1  modules. 

Surfaces 

Surfaces  are  modeled  in  finite-element  fashion, 
currently  as  contiguous  triangular  elements  joined  at 
the  vertices  (nodes).  PIMS  does  not  generate  3-D 
geometries  or  surfaces;  instead,  it  imports  them  from 
other  software. 

Users  create  custom  geometries  using  almost  any 
mainstream  commercial  3-D  solid  modeling  package. 
Then,  they  use  finite  element  analysis  software  to 
mesh  the  surface  of  their  geometry  as  if  they  were 
going  to  perform  a  structural  analysis  using  thin  shells. 
The  user  then  saves  the  meshed  surface  file  in  ANSYS 
format,  which  is  readable  by  PIMS.  ANSYS  finite 
element  format  was  chosen  because  it  is  widely 
supported  by  finite  element  packages. 

This  concept  of  separating  the  surface  geometry 
definition  from  the  plasma  calculation  has  proven  very 


Fig.  2.  Sputter  yields  for  bombardment  by  singly 
ionized  xenon  at  300eV. 

successful.  It  greatly  reduced  development  time  and 
cost  by  eliminating  the  need  for  a  separate  surface 
definition  module.  It  allows  users  to  choose  which 
software  to  use  in  defining  geometries.  And,  users  can 
import  into  PIMS  geometries  that  have  already  been 
defined  for  other  reasons  (structural,  thermal,  etc.). 

Material  Properties 

The  user  constructs  a  database  of  materials  that  is  read 
by  PIMS.  The  database  contains  material  names, 
material  reference  numbers,  and  molecular  weights 
and  charges  (in  the  case  of  ions).  Materials  in  the 
database  are  connected  to  the  surface  geometry  by  the 
material  reference  number.  Users  mark  surface 
materials  during  geometry/surface  definition  using 
their  finite  element  software.  They  simply  set  the 
elastic  modulus  of  the  surface  component  to  be  equal 
to  the  material  reference  number.  This  value  appears 
in  the  ANSYS  file,  where  PIMS  can  read  it. 

The  user  also  provides  a  second  database,  a  materials 
interaction  database.  This  database  contains  the 
sputter  yield  coefficients  and  sticking  coefficients  of 
one  material  interacting  with  the  other.  For  example, 


one  important  interaction  may  be  between  Xe+  and 
Kapton. 

The  following  equation  has  been  used  to  generate  the 
sputter  yield  for  each  species  as  a  function  of  the  ion 
energy  and  incidence  angle.6 

Y{E,e)  =  {a  +  bE)(\.0-0J2e  +  \\.12e2 -lAW3 -2.5ieA) 

Above,  E  is  the  particle  energy,  and  9  is  the  incidence 
angle  (off-normal).  For  materials  considered  here, 
Table  1  gives  the  coefficients  used,  and  Fig.  2  plots 
sputter  yield  for  E  =300eV. 

Sources 

Sources  are  modeled  as  having  a  specific  velocity 
distribution  that  is  constant  over  individual  surface 
elements.  A  collection  of  commands  allows  the  user 
to  either  specify  one  of  a  set  of  source  types  (mono- 
energetic,  half-Maxwellian,  etc.)  or  read  in  a  file 
containing  a  custom  discretized  velocity  distribution 
function. 

Source  elements  are  identified  with  a  source  reference 
number  during  geometry/surface  definitions,  much  like 
the  materials  are  identified. 

Although  this  method  is  extremely  descriptive  and 
general,  Level  1  modules  may  treat  this  information  in 
various  ways.  For  instance,  a  Level  1  module  could  be 
written  to  treat  the  source  element  a  single  point 
source  for  ray  tracing  purposes.  Alternately,  particle 
methods  could  sample  from  the  velocity  distribution 
and  introduce  particles  over  the  full  element  surface. 
Therefore,  this  choice  of  source  definition  methods 
gives  PIMS  the  greatest  flexibility. 

Plasma  Simulation 

Currently,  two  Level  1  modules  have  been  written. 
The  first,  PRESCRIBED JPLUME,  allows  the  user  to 
import  a  previously  calculated  or  measured  plume 
field.  This  plume  is  superimposed  over  the  user’s 
surface  geometry.  Plasma  densities,  fluxes,  and 
sputter  rates  are  then  calculated  at  each  surface  node. 

The  second  module,  RAY,  uses  ray  tracing  to  calculate 
the  flux  from  all  sources  onto  all  surface  nodes.  Once 
again,  density,  flux,  and  sputter  rate  are  calculated. 


#  pirns . in 

# 

#  Load  a  GEO  satellite  geometry,  add 

#  a  3kW  HET  source ,  calculate  the  flux 

#  and  sputtering  using  ray  tracing, 

#  and  save  the  results  in  Tecplot  format. 

# 

surf ace__load  ANSYS  GEO_Sat.ANS 
source  HET  27  het_J3kw.dat 
ray 

surface_save  TECPLOT  GEO_Sat .  dat _ _ 


Fig.  3.  Sample  PIMS  command  file. 


Future  modules  will  incorporate  statistical  kinetic 
methods  for  plasma  calculation  such  as  PIC  and 
DSMC.  Plans  also  include  development  of  kinetic 
algorithms  for  use  on  unstructured  meshes,  adaptive 
meshes,  and  domain  decomposition.  Primarily,  these 
techniques  will  be  incorporated  to  add  flexibility  to  the 
simulation.  For  instance,  domain  decomposition  will 
allow  the  domain  to  be  broken  into  smaller  sub- 
domains,  each  potentially  having  different  algorithms, 
depending  on  local  parameters  as  the  Debye  length  or 
mean  free  path. 

User  Interface 

The  user  enters  commands  via  a  PIMS  input  file.  The 
commands  are  executed  sequentially  as  they  appear  in 
the  input  file.  Each  command  may  have  some  number 
of  parameters  separated  by  spaces  or  commas.  A 
sample  input  file  is  shown  in  Fig.  3. 

Typical  run  times  for  low-fidelity  cases  (using 
PRESCRIBED  PLUME  or  RAY)  take  approximately 
20  seconds  on  a  2  GHz  Intel  Pentium  4  workstation. 
Once  more  detailed  physics  are  incorporated,  with 
Level  1  modules  incorporating  such  algorithms  as 
PIC-DSMC,  run  times  are  expected  to  be  between  20 
minutes  and  20  hours,  depending  on  the  level  of 
fidelity  and  on  the  initial  conditions. 

Results  and  Discussion 

Initial  test  PIMS,  runs  were  executed  for  two  cases:  a) 
a  fictitious  geosynchronous  satellite  with  an  HET 


firing  in  the  north  direction  (as  if  for  stationkeeping), 
and  b)  an  HET  firing  inside  a  vacuum  chamber  (as  if 
during  a  flight  readiness  test).  In  both  cases,  the  Level 
1  module,  PRESCRIBED_PLUME  was  used  to 
incorporate  a  previously  calculated  plume  expansion 
model  onto  the  surface  geometry.  The  plume 
expansion  model  used  here  was  calculated  for  a  Busek 
200-Watt  HET7  by  SAIC  using  the  GILBERT5 
toolbox. 

Results  from  the  first  case  are  shown  in  Fig.  4  through 
Fig.  7.  Fig.  4  shows  the  geometry  of  the 
geosynchronous  satellite  model,  with  surfaces  broken 
down  into  triangular  elements.  Two  commercial 
packages,  SolidWorks  and  COSMOSWorks,  were 
used  to  generate  these  surface  geometries.  Final 
plotting  was  performed  by  another  commercially 
available  package,  Tecplot.  The  colors  of  the  mesh 
lines  indicate  the  type  of  material.  Fig.  5  shows  a  slice 
through  the  200- Watt  HET  plume  superimposed  on  the 
satellite  model.  Plasma  density  is  highest  near  the 
HET  exhaust,  and  drops  off  rapidly  as  the  plume 
expands  upward  toward  the  solar  arrays.  Fig.  6  shows 
the  PIMS  calculation  of  ion  flux  at  the  surface  of  the 
satellite.  Finally,  PIMS  calculates  the  rate  of  surface 
sputtering,  given  in  Fig.  7. 

The  sputtering  rate  peaks  near  the  solar  array  comer. 
This  illustrates  a  real  problem  with  electric  propulsion 
on  geosynchronous  satellites.  For  north- south 
stationkeeping,  the  ideal  firing  direction  (from  a  thrust 
efficiency  standpoint)  is  directly  north.  However, 
PIMS  shows  us  that  long-term  firing  of  the  HET  over 
the  lifetime  of  a  satellite  (~7000  hours)  in  this 
configuration  may  remove  2.5  mm  from  the  surface  of 
the  solar  array  at  the  corner.  In  reality,  the  solar  array 
will  be  rotating  to  track  the  sun,  and  will  not  always 
have  a  comer  directly  in  the  HET  plume.  So  2.5  mm 
can  be  considered  a  worst  case.  Other  ways  of 
reducing  the  sputtering  are  to  angle  the  HET  plume 
away  from  due  north. 

Results  from  the  second  case  are  shown  in  Fig.  8 
through  Fig.  11.  Fig.  8  shows  the  geometry  of  the 
interior  of  a  vacuum  chamber,  with  an  HET  in  the 
center,  and  several  plasma  diagnostic  instruments 
arrayed  nearby.  Once  again  the  200- Watt  HET  plume 
is  incorporated,  and  fluxes  and  sputtering  rates  are 
calculated  by  PIMS. 


Referring  to  Fig.  10,  a  spherical  plasma  probe  can  be 
seen  protruding  from  the  instrumentation  panel.  Ion 
flux  to  this  probe  can  be  taken  directly  from  the  figure. 
The  peak,  on  the  HET-side  of  the  probe,  is 
approximately  lxlO18  mV1. 

Fig.  1 1  shows  an  interesting  pattern  of  sputtering  on 
the  electron  probe.  Maximum  sputter  rate  appears  as  a 
ring  around  the  point  on  the  probe  nearest  the  HET. 
This  is  due  to  the  dependency  of  sputter  yield  on 
incidence  angle.  Sputter  rate  at  normal  incidence  is 
typically  lower  than  that  at  grazing  angles,  as  can  be 
seen  in  the  discussion  of  sputter  models  above. 

Conclusions 

Although  still  in  an  early  stage  of  development,  The 
Plasma  Interaction  Modeling  System  (PIMS)  now  can 
help  predict  ion  flux  and  sputtering  of  surface 
materials  both  onboard  spacecraft  and  in  laboratory 
test  facilities.  PIMS’  modular  architecture  is  allowing 
rapid  expansion  of  its  capabilities,  and  giving  users 
flexibility  to  design  their  own  geometries  and  choose 
their  preferential  plasma  simulation  method. 

Additional  work  for  the  future  includes  expansion  of 
the  source  module,  incorporation  of  surface  re¬ 
deposition,  and  construction  of  new  Level  1  modules 
that  can  self-consistently  compute  plasma  expansion 
and  interaction  with  surfaces. 
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Fig.  6.  Flux  of  xenon  ions  to  the  satellite. 


Fig.  7.  Surface  sputtering  rate. 


Fig.  10.  Flux  of  xenon  ions  to  the  surfaces. 


Fig.  11.  Surface  sputtering  rate. 


